arXiv:1508.06968vl [cond-mat.stat-mech] 27 Aug 2015 


Correlation functions in conformal invariant stochastic 

processes 

Francisco C. Alcaraz^ and Vladimir Rittenberg^ 


^Instituto de Fisica de Sdo Carlos, Universidade de Sdo Paulo, Caixa Postal 369, 

13560-590, Sdo Carlos, SP, Brazil 

Physikalisches Institut, Universitdt Bonn, Nussallee 12, 53115 Bonn, Germany 

August 28, 2015 


Abstract 

We consider the problem of correlation functions in the stationary states of one¬ 
dimensional stochastic models having conformal invariance. If one considers the space 
dependence of the correlators, the novel aspect is that although one considers systems 
with periodic boundary conditions, the observables are described by boundary oper¬ 
ators. From our experience with equilibrium problems one would have expected bulk 
operators. Boundary operators have correlators having critical exponents being half 
of those of bulk operators. If one studies the space-time dependence of the two-point 
function, one has to consider one boundary and one bulk operators. The Raise and Peel 
model has conformal invariance as can be shown in the spin 1/2 basis of the Hamilto¬ 
nian which gives the time evolution of the system. This is an XXZ quantum chain with 
twisted boundary condition and local interactions. This Hamiltonian is integrable and 
the spectrum is known in the finite-size scaling limit. In the stochastic base in which 
the process is defined, the Hamiltonian is not local anymore. The mapping into an SOS 
model, helps to define new local operators. As a byproduct some new properties of the 
SOS model are conjectured. The predictions of conformal invariance are discussed in 
the new framework and compared with Monte Carlo simulations. 


1 Introduction 

Many-body local one-dimensional stochastic processes have received a lot of attention in 
both experimental [1] and especially theoretical research. A special class of models are 
those who have conformal invariance as a symmetry. They have a dynamic critical exponent 
z = 1. Conformal invariance gives a lot of information about the expected behavior of 
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the system in the hnite-size scaling limit (large sizes and times) describing the approach 
to the stationary state) |2]. Another important prediction of conformal invariance is about 
the expected behavior of the correlation functions at criticality. One takes the size of the 
system to inhnity hrst and look at the large distance behavior of the correlators next. For 
equilibrium systems the problem is well understood but this is not the case for stochastic 
processes and it is the aim of this paper to clarify it. 

We will explain why the generalization to non-equilibrium processes of correlators is 
not a trivial matter. When studying one-dimensional equilibrium quantum systems with 
periodic boundary conditions, one takes average quantities in the ground-state of the system: 
( 0 |F(a;i, ti)F(a; 2 , t 2 )| 0 ) where F is a local operator (in general a primary field), | 0 ) = ( 0 |’'' 
is the ground-state wave function of the Hamiltonian which gives the time evolution of the 
system. For simplicity, one can see the Hamiltonian as describing a quantum chain. The 
space-time behavior of the correlator for a periodic system is 

{0\F{x^,h)F{x2,t2)\0) = ( 1 . 1 ) 

where = (x 2 — Xi)^ + (^2 — H)^ and Xb is a critical exponent {xb = A + A, where A and 
A are the scaling dimensions). The scaling dimensions are known from the representations 
of the Virasoro algebra which describe the system in the continuum space and time in the 
hnite-size scaling limit. The space and time dehne a complex plane or an inhnite cylinder of 
a very large radius. For an open system one considers a half-plane respectively a strip. 

The situation is very different in the case of stochastic processes. We remind the reader 
that in the continuum time description, the Hamiltonian with matrix elements Ha,b (o, b = 
1,2,, N) is special: the non-diagonal elements Ha,b are non-negative (they are the rates 
for the transitions b a) and the diagonal ones are hxed by the non-diagonal ones: Ha,a = 
— J2b,b^a Ha,b- A matrix with this property is called an intensity matrix and we will call 
the basis in which H is an intensity matrix, a stochastic basis. We stress that a change of 
basis through a similarity transformation changes in general, the matrix from an intensity 
matrix to an usual matrix. In the stochastic basis, the bra vector (0|^ = (1,1,1,..., 1| has 
the property {0\^H = 0, for any value of N. The time evolution of the system is given by 
the master equation: 

= ( 1 . 2 ) 

Pa(t) gives the probability to hnd the state a at the time t. For the stationary state, the 
probabilities give a ket-vector | 0 ) = {P^, P^, ..., P^) which is an eigenvector: F| 0 ) = 0 . 
Notice that {0|^ 7 ^ |0)+. 

To clarify the picture, it is useful to have the following scenario in mind. Assume that 
one has a non-stochastic basis |a) (a = 1, 2,..., N) in which H is Hermitian and describes 
an integrable quantum chain with local interactions and periodic boundary conditions. We 
also assume that the spectrum of H is known and that in the hnite-size scaling limit is given 
by Virasoro representations with a central charge c = 0 (keep in mind that in a stochastic 
process the ground-state energy is zero for any system size). Any correlator of local operators 
can be computed in the standard way and one obtains expressions like (II.Ih . If now we make 
the change of basis to the stochastic one, the spectrum stays unchanged, but the action of 
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the Hamiltonian in the stochastic basis is not local anymore but stays periodic. We will 
come back to this important point later in the text. The dehnition of a correlator in the 
stationary state is: 

(0|^G(a:i, ti)G(x2, ^ 2 ) 10 ), (1-3) 

where G{x,t) is a diagonal operator in the stochastic basis. Since < 0|^ is a trivial vector, 
the expression fll.Sp looks more like a form factor than a correlator. 

To simplify the problem, let us look in the stationary state at space correlators only. The 
way to solve the problem can be found in a seminal paper by Jacobsen and Salem [3]. Their 
interest was to compute some properties of the ground-state wave function of the XXZ chain 
in the context of information theory. They considered the complex half-plane with x on the 
horizontal and y on the vertical axis [y > 0). On the x-axis one considers two boundary 
helds |1] G(a;i) and G{x 2 ), and assumes that one has Neumann boundary conditions on the 
x-axis (this corresponds to the bra vector (0|^ in the quantum chain). One makes a conformal 
transformation of the half-plane to a half-cylinder (a cylinder with a cut perpendicular to 
the cylinder axis). The x-axis is transformed in the circle which cuts the cylinder with a 
coordinate u and the |/-axis is transformed in the t (time) half axis along the cylinder. The 
time t evolution of the system is given by the Hamiltonian in the stochastic basis. One 
assumes that the initial state of the stochastic system is taken at t = —00 and that the 
system ends up in the stationary state at f = 0 (the circle at the cut). The boundary 
operators G become the observables in the stochastic process and their correlator can be 
computed using the conformal held results for boundary operators in the half-plane. Since 
the correlation functions of boundary operators have the critical exponents given by only 
one Virasoro algebra and not two, like for bulk operators, it implies, roughly speaking that 
the critical exponents seen in stationary states have half the values of those of equilibrium 
states. 

In Sec. 2 we extend this observation to the non-periodic open system and to time corre¬ 
lators. To compute the latter, one needs to consider boundary and bulk operators. 

We come now to the applications of conformal invariance to stochastic models. Unfor¬ 
tunately very few are known and they are all versions of the Raise and Peel model mu 
that we discuss in this paper. The model is reviewed in Sec. 3. One starts with the periodic 
Temperley-Lieb algebra with the parameter q = exp(i7r/3) and write the Hamiltonian as a 
sum of its generators. The algebra has several representations which are considered in the 
text. In the spin 1/2 representation one obtains the XXZ quantum chain with an anisotropy 
A = 1/2 and twisted boundary conditions. The Hamiltonian is Hermitian and its spectrum 
is massless and given by known representations of the Virasoro algebra with a central charge 
c = 0, hence conformal invariant. The stochastic basis is obtained using other representa¬ 
tions of the algebra. We use three stochastic base. The hrst one is the link patterns one, the 
second is the Dyck paths, the third is the particle-vacancy basis. They are equivalent. For 
example, the particle-vacancy representation is obtained taking the slopes of the Dyck paths. 
A positive slope is a particle, a negative slope is a vacancy. In these representations, the 
Hamiltonian acts non locally and it is not obvious to hnd observables which can be identihed 
with boundary or bulk operators. Although it is relatively easy to hnd observables having 
the proper space dependence which can be compared with the predictions of conformal in¬ 
variance it might be very difficult, if not impossible, to get operators with the proper time 
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dependence. In other words, it is much easier to identify local boundary operators and much 
harder to identify bulk operators. To take an example, the current density was computed 
and a conjecture for any lattice size was given in [7]. For a given conhguration, the current 
density gets contributions from the action of the generators of the Temperley- Lieb algebra 
on all the sites of the system, this makes probably, impossible to look for a proper time 
dependence which describes a local behavior. The reason is the following one. In the time 
interval At/L which is the time scale for the continuum limit (sequential updating in Monte 
Carlo simulations) only one generator hits the conhguration and not all of them. Another 
problem is related to the computation of the correlators using Monte Carlo simulations in 
the domain t/L « 1 (L is the system size). For the accessible values of L, the times are so 
short that one ends up with results depending on the initial conditions. 

In the case of the Raise and Peel model, one could hnd local operators and their space 
dependence, looking at the SOS (Coulomb gas) model picture of boundary vertex operators. 
This is the topic of Sec. 4. Based on the work of |1], we could identify one operator (called 
the NCA for non-crossing arches) and the contact point operator. This identihcation was 
possible using the generating function for counting the total number of arches crossing the 
ends of a segment. As a by-product of our research, we present a conjecture for the generating 
function giving the probabilities to have n crossings at one end and m at the other end of the 
segment. We also show, using the particle-vacancy stochastic basis that the particle density 
operator is also a ’’good” boundary operator. 

Once the boundary operators have been identihed, in Sec. 5 we compare the predictions 
of conformal invariance with the data obtained from Monte Carlo simulations. We also show 
that unlike the NCA two-point correlation function, the two-contact point correlation has 
a simple analytic expression but vanishes in the large L limit. The problems which appear 
when dealing with non-locality are outlined. 

In Sec. 6 we present our conclusions. 


2 Conformal invariance and stochastic processes 

In this section we present the predictions of conformal held theory for the one and two-point 
correlation functions seen in the stationary state of a stochastic process. These predictions 
will be compared with various correlators in Secs. 4 and 5. 

Consider a one-dimensional system with N states (conhgurations) with probabilities Pa{t) 
{a = 1,... ,N). The evolution of the system in continuum time, is given by the master 
equation fll.2p . 

The Hamiltonian H is an intensity matrix (dehned in the Introduction). One gives the 
probabilities at the initial time, say t = —oo, and the system evolves to a stationary state 
which is reached at the time t = 0. One is interested to know the correlation functions 
in the stationary state. We assume that the N states can be seen as conhgurations of a 
one-dimensional lattice with L sites and that the spectrum of H was computed and found 
that the hnite-size scaling spectrum could be organized in representations of the Virasoro 
algebra and that hence one has conformal invariance. The central charge of the Virasoro 
algebra is c = 0 since in a stochastic process the ground-state energy is equal to zero for any 
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Figure 1: Periodic boundary conditions, a) The observables in the stationary state of a 
stochastic process G{ui) and G{u 2 ) are attached to the rim of the half-cylinder. They are 
used when computing the space two-point function. B{u,t) is an observable on the surface 
of the half-cylinder. It is used when computing space and time correlators, b) G{xi) and 
G{x 2 ) are boundary helds in the half-plane, c) G(0) is a boundary held at x = 0, 5(0,?/) is 
a bulk held at x = 0 and y. The Hamiltonian H at time t (negative) acts on the circle of 
length L 

number of sites. We would like to know what are the consequences of conformal invariance 
on the correlation functions fll.31) . The latter can be obtained, as discussed in Secs. 4 and 5, 
through computer simulations in which one discretize the time by taking Monte Carlo steps 
At/L or by using held theory (taking the lattice space to zero). 

In order to hnd the consequences of conformal invariance on the correlation functions, 
we take the continuum space and time limits and consider separately the cases of periodic 
and open boundary conditions. 

a) Periodic boundary conditions. 

To start with, we consider a half-cylinder with its half-axis being the time {t < 0) 
coordinate (see Fig. la). The cylinder is cut at f = 0 which corresponds to the time when 
the system reaches the stationary state. The horizontal circles on which H is dehned, have 
a perimeter L. We introduce a second coordinate u on the circles (0 < u < L). The 
systems evolves from the far away bottom of the half-cylinder (the initial state) to its top 
(the stationary state). The observables of the stochastic process in the stationary state 
{t = 0) used in the calculation of the space dependent two-point functions will be denoted 
by G{ui) and G{u 2 ) (they are on the top of the half-cylinder). If one studies time and space 
two-point correlations, one observable G{ui) is on the top of the half-cylinder and the second 
one on its surface. 

In order to get the space correlation functions on the half-cylinder, one uses the fact that 
they are known in the half-plane and that one can make a conformal transformation from 
the half-plane to the half cylinder. We denote hy z = x + iy {y > 0) and w = t + iu (t <0) 
the complex coordinates, and take Neumann boundary conditions aX y = 0. This choice 
corresponds to (0|® = (1,1,...,!) in the stochastic process. The correlation function of two 
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boundary fields G{xi) and G{x 2 ) (they are situated a.t y = 0 and shown in Fig. lb) is known 
from conformal field theory [S]: 

(0|G(a;i)G(a;2)|0) = jx2 - (2.1) 


Here A is the scalling dimension of the boundary field G. The Kac table provides the 
following possible values for a c = 0 theory: 


A 


(3r - 2sf - 1 
24 


( 2 . 2 ) 


where r and s are nonnegative integers (A = 0,1/8,1/3, 5/8,1, 35/24,...). 
The transformation 



/I + iz\ 
[l-izj 


(2.3) 


takes the half-plane to the half-cylinder and the two-point correlation function on the top of 
the half-cylinder follows: 


(0|-G(«,)G(«2)|0) = (0|G(xOG(x.)|0). 

Taking L —)■ cxo with ui and U 2 fixed, as well |m 2 — Ui\ large, one obtains: 


(0|*G(Mi)G(n2)|0) = \u2 - ui 


(2.4) 


(2.5) 


This is an interesting result: although the stochastic process takes place in a system with 
periodic boundary conditions, one doesn’t have left and right movers. One has only one Vira- 
soro algebra and not two as it is the case for periodic boundary conditions in the equilibrium 
problem. 

The one-point function vanishes for large L. The large L behavior is given by the scaling 
dimension of the boundary operator. Assuming that on the half-line. 


{0rG(a:)|0) = G, (2.6) 

where G is a constant and making the mapping fl2.3p one obtains, for large L, the one-point 
function on the rim of the cylinder: 


Gl{u) ~ L-^. (2.7) 

In Sec. 4 and 5 we are going to make repeated use of Eq. fl2.7p . 

If one is interested in time correlation functions (0|^G(0)i?(0,t)|0) (for simplicity both 
observables have u = 0), one looks again at the half-plane (see Fig. Ic) taking x = 0 (i. e. 
G(x = 0)) and B{x = 0,?/). Like before, G is a boundary field and i? is a bulk field. We 
have to keep in mind that the limit ?/ —)■ 0 is singular (the bulk field B doesn’t go smoothly 
to a boundary field). One has |1]: 

(0|G(0)i?(0,|/)|0) = 2/-"--^, (2.8) 
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Figure 2: Open system, a) the strip w = t + iu on which the stochastic process takes place. 
The stationary state is at t = 0. The boundary operator G{u) is attached on the top of the 
strip, b) The half-plane z = x + iy in which the conformal held theory is used. C'(—1) and 
0 (1) are boundary changing helds applied at x = — 1 and x = 1. 

where is the conformal dimension of the bulk held B. We perform again the conformal 
transformation from the half-plane to the half-cylinder and get: 

(0|"G(m = 0)5(m = 0,t)|0) = (2.9) 

In order to derive fl2.9l) we have taken L ^ oo with t hxed. The expression fl2.8p generalizes 
in a simple way when G and B have diherent space coordinates, say Ui and U 2 - 

b) Open boundaries. 

If the stochastic process has open boundaries, the problem is slightly more complicated. 
The system is not dehned on a half-cylinder but on a half-strip of width L (see Fig. [2]). 
We consider the case of only one observable in the stationary state and study the hnite-size 
scaling ehects on the one-point function Gl{u) {u = 0,... ,L). The two-point correlation 
function implies the knowledge of model dependent correlations in the half-plane. 

In order to get the expression of the one-point function, one consider the half-plane with 
the 1 / = 0 axis divided by the points A {x = —1) and O (x = 1). Between A and O ( 

— 1 < X < 1) one takes Neumann boundary conditions. For x < —1 and x > 1, we don’t 
need to specify the boundary conditions except that they are the same. We attach to the 
points A and O, boundary changing helds |lj 0(—1) and 0(1). In order to hnd 0(x), we 
consider the vacuum expectation of the 3-point function given by conformal hied theory: 

{0\<i>a{Xi)<i,t(x2)<l>c(x3)\0) >= _ j,3)A,+A,-A. )A.+A.-A, ■ 

(2.10) 

where Aq, Af, and Ac are the scaling dimensions of 0a, 06 and 0c. Taking Aq = Ac, Af, = A 
and using (12.101) one gets: 

(0|0 (-l)0(x)0(l)|0)~(l-x^)-^. (2.11) 

We use the Schwarz-Christohel mapping which takes the half-plane to the strip: 2 ; = 

— cosh(7rtc/L) {z = x + iy,w = t + iu) and obtain 

Gi(«)=(^) {0|C(-l)G(i)C(l)|0)~ (Ism(Wi)) (2.12) 
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This result is interesting since there is a theorem [S] which says that in equilibrium, the 
one-point function vanishes unless one has a left-hand mover with a scaling dimension equal 
to the one of the right-hand mover therefore one has: 

L 

— sin(7rM/L)) . (2.13) 

71 } 

Notice the exponent A in (I2.12p as compared to 2A in fl2.13p . 

For not symmetric boundary conditions, similar to the equilibrium problem [9], Eq. (I2T2]) 
generalizes as follows: 

/ T \ 

Gl(w) ~ E(cos( 7™/L)) f — sin(7TO/L) j . (2-14) 

The time dependence of the one-point function can also be easily obtained but for the 
application we have in mind (I2.12p is enough. 

We have to stress that these results were obtained assuming that in the stochastic process 
one uses observables which are local. As we are going to see in the next sections, to hnd 
local observables is not an easy task. 

In Sec. 4 and 5 we compare the expressions fl2.5p . fl2.7p . fl2.12p and fl2.14p with the 
correlation functions seen in the Raise and Peel model. This model is described in the next 
section. 



3 The Raise and Peel model 

There are plenty of papers [5] describing this model, written after the discovery by Razumov 
and Stroganov HU, that its stationary state has magic combinatorial properties. We give 
here a short survey of the model which should allow the reader to follow the results presented 
in Secs. 4 and 5. 

One considers a Hamiltonian given by the expression: 

M 

» = ( 3 ' 1 ) 

i=l 

where M = L — 1, for the open system and M = L for the periodic one. ej(i = 1,..., M) 
are the generators of the Temperley-Lieb algebra for the open system and of the periodic 
Temperley Lieb algebra for the periodic one. They verify the relations 

e? = {q + l/q)ei (3.2) 

e*efc = CfeCi, \i-k\>l. (3.4) 

In the case of the periodic algebra, which is inhnite dimensional, another relation is 
needed in order to get a hnite dimensional quotient m- We skip it here. The algebras have 
a representation in terms of Pauli matrices and one obtains a Hamiltonian given by an XXZ 
























Figure 3: The three stochastic base, Dyck paths, particle-vacancy and link patterns base 
for L = A. There are 6 possible prohles. In the particle-vacancy representation the particles 
are full dots and the vacancies empty dots. The corresponding link patterns representation 
are also shown. 

quantum chain with L sites having the Ug{sl{2)) symmetry for the open system. In this case 
the Hamiltonian is not Hermitian. In the periodic case, one gets the XXZ chain with twisted 
boundary conditions where the twist depends on q. H is Hermitian in this case. The spectra 
of the Hamiltonians are known and for |g| = 1, they are gapless. Taking q = exp(i7r/3) the 
ground-state eigenvalue is zero for any value of L, and the other energy levels are positive 
and therefore one can look for stochastic base. 

We consider the case L even only. There are three equivalent stochastic base, each one 
useful for different purposes. The simplest one is given by Dyck paths. They are dehned in 
terms of heights {hi} obeying the restricted solid-on-solid rules: 

hi+i - hi = ±1, i = 0,1,. .. ,L, (3.5) 

with ho = hi = 0 for the open system and hi = hL+i for the periodic one. In Fig. [3] we 
show the six conhgurations in the periodic case L = 4. There are L\/{{L/2 + l)([L/2]!)^} 
conhgurations for the open system and L!/[(L/2)!]^ for the periodic case. 

In the stochastic basis of Dyck paths the time evolution of the system can be visualized 
in the following way. A Dyck path can be seen as separating droplets of tilted tiles from a 
rarehed gas of tilted tiles (see Figure |4]). The Dyck path changes its shape as a result of hits 
from the tiles of the gas according to the following rules: 

In a time interval At, with probability pi = At / L a. given site i, with height hi is reached 
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1 2 18 

Figure 4: Example of a configuration with four peaks and three contact points in the 
periodic Raise ans Peel model with L = 18 sites. Depending on the position where the tilted 
tiles reach the interface, several distinct processes occur (see the text). 


by a tilted tile of the gas phase (see Fig. 0]). Several movements are possible depending on 
the local slop s* = (hj+i — hi)/2 at the site: 

i) Si = Q and hi > hi-i. 

The tile hits a local peak (case a in Fig. 0)) and it is reflected leaving the prohle unchanged. 

ii) Si = 0 and hi < hi_i. 

The tile hits a local minimum (case b in Fig. 0]). The tile is absorved {hi hi + 2). 

Hi) Si = 1. 

The tile is reflected after triggering the desorption of a layer of tiles from the segment 
hj > hi = hi+b, j = i + 1,... ,i + b — 1, i. e., hj ^ hj — 2 ioT j = i + 1,... ,i + b — 1 (case c 
of Fig. 0]). The desorbed layer contains 6 — 1 tiles (always an odd number). 

iv) Si = —1. 

The tile is reflected after triggering the desorption of a layer of tiles from the segment 
hj > hi = hi+b, j = i — b + 1,... ,i — 1, i. e., hj ^ hj — 2 for j = i — b + 1,... ,i — 1 (case d 
of Fig. 0]). 

These evolution rules take place in discrete time by choosing, for example. At = 1 as 
usual in Monte Carlo simulations or in continuum time by taking At = 1/L —)■ 0. 

There is a one-to-one correspondence among the Dyck path conhgurations {hi} and the 
conhgurations {rii} of particles {ui = 1) and vacancies (rij = 0) attached to the links of the 
L sites chain, namely, 

rii = ' (* = 1>2,... ,L), hL+i = hi. (3.6) 

For illustration the conhgurations for the L = 4 sites system are also shown in Fig. 01 
This correspondence imply that the number of particles and vacancies are equal to L/2. 
The time-evolution rules of this half-hlling chain of particles are directly obtained from the 
corresponding rules for the heights we dehned (see Cl)- In this particle language we have 
a nonlocal asymmetric exclusion model where the particles either jump to the next-nearest 
neighbor position (corresponding to adsorption) at its left, or jump nonlocally to a site 
position on its right (corresponding to desorption). 

We have just described the second stochastic basis, the particle-vacancy one. This repre¬ 
sentation as well as the Dyck path one correspond to the = 0 sector of the XXZ quantum 
chain. It can be enlarged to the full 2^ dimensional representation as shown in [7]. A third 
stochastic basis is the link patterns basis. It is obtained from the Dyck paths basis in the 
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following way. For each height hi {i = 1, ...,L) one draws hi noncrossing arches which end 
up at different sites. When all sites i are considered, arches end up at all sites. In Fig. 4 we 
illustrate the correspondence between the two base for L = 4. The dynamics of the system 
in the link patterns basis follows the one in the Dyck path basis. 


4 Connection with the loop gas and the SOS models 

The stationary state of the Raise and Peel model with periodic boundary conditions in the 
link patterns basis can be understood in terms of a Q = 1 loop gas model |1] on a half¬ 
cylinder (see Figs. [I]and|5]). We have in mind the discussion in Sec. 2 about the applications 
of conformal invariance. The rim of the cylinder corresponds to the stationary state. We 
give here a short summary of this connection (see for more details). 

One takes a lattice on the half-cylinder and consider only the conhgurations of self and 
mutually avoiding fully packed loops. The conhgurations of the loop model are either closed 
loops (they don’t reach the rim of the cylinder at t = 0) shown in black, or arches which 
start and end on the rim, colored in Fig. [5l All arches have a fugacity equal to 1. The 
probability to have a certain conhguration in the stationary state of the stochastic model is 
proportional to the partition function on the half-cylinder with the given conhguration as a 
boundary condition. This observation is relevant since one can map the loop model into a 
solid-on-solid (SOS) model on the lattice. In the continuum space and time limit, the latter 
model is described by a free bosonic conformal held theory on the rim of the half-cylinder 
[1]. The model is parameterized by the parameters eo equal to 1/3 and the coupling g = 2/3 
in our case. If we denote by G{u) the boundary bosonic held (with Neumann boundary 
conditions on the boundary) the correlator is 


< G{u)G{u') >= —3/2 In \u — u'\^. 

(4.1) 

We dehne the vertex operators V±{u) = exp(±i(ei -f eQ/2)G{u)). The 
parameter ei will be clear soon. The correlator of two vertex operators is 

use of the new 

< V/_(m)I/_(m') >~ \u — u' 

(4.2) 

with 

A = (36e? - l)/24. 

We can now use Eos. 112.61). 112.711. 114.311 and eet for laree values of L\ 

(4.3) 

< 14 >= A(ei)L-^. 

(4.4) 

A{ei) is a function which has to be determined. Consider now the generating function 

M(w) = 

(4.5) 


{hi} 


where P{hi} is the probability to have a height hi (Dyck path representation) at the site 
i. The site independence of M{w) is a consequence of the translational invariance of the 
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stochastic model. This corresponds to hi arches crossing the link i between the sites i and 
i + 1 (see Fig. |H]) 

It is important to note that following 013] we can get the identity 

< 14 >= M{w), (4.6) 

where 

tc = 2cos(7rei)/\/3. (4.7) 

In the next section we will discuss the implications of the relation 

M{w) = A{w)L-^^^^ (4.8) 


for the Raise and Peel model. Although the relation fl4.8p might have been derived somewhere 
else, we have checked it using computer simulations. In Fig. [7] we show the values of the 
generating function M{w) as a function of the lattice size L for 3 values of w. 


wi = ^, ^ = \ = 0.125; W2 = \ % A = 0.05208333...; 

\/3 8 V 3 

\/3 

W3 = ^, A = 0.037720222... . 


The hts are done using Eq. fl4.8p and they are excellent. From the hts we obtain 
>1(1/73) = 0.850; A{V2/3) = 0.936; >1(73/2) = 0.954. 


(4.9) 


(4.10) 


The value 4(1) = 1 is a consequence of the normalization of the probability function 
in Eq. fl4.5p . For completeness we also give the value 

4(0) = 2^^^^ = 0.424... (4.11) 

3y TT 

This value will be derived in the next section. 

We consider now the two-point function fl4.2p and look at the lattice realization of it. We 
take a segment {i,j) on the rim of the half-cylinder and consider the generating function: 

M'{u;w) = ^ n;"“+"^P(4,4), (4.12) 

hi^hj 

where P(4, hj) is the probability to have the heights hi and hj, and u = j — i. If one looks 
at the computer simulations one discovers that the function M'{u]w) is not independent on 
L for large lattice sizes. To hnd the proper dehnition of the generating function, it is better 
to use the arches language. If one considers the segment (see Fig. [6]), there are two sorts of 
arches. The hrst one are those which start at sites smaller than i and end up at sites larger 
than j or starts and ends at the sites between i and j. The second one are the arches which 
start inside and end outside the segment between i and j. We call them non-crossing arches 
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Figure 5: A configuration of self and mutually avoiding fully packed loop model with 
L = 12 sites in the periodic spacial direction (horizontal). The time direction is the vertical 
direction. The closed arches in the bulk (black loops) as well the closed arches (colored 
arches) that starts and close at the surface (time t) have a Boltzmann weight 1. 



Figure 6: A configuration with the segment hi = 2, hj = 3. One arch enters the 

segment at i and leaves it at j. 

(NCA). If we disregard the hrst kind of arches, in the calculation in the new generating 
function M{w) one hnds, in agreement with fl4.2p 

M{u- w) = (4.13) 

with A given by (I4.3p . This implies B[l) = 1, again from the normalization of the probability 
function. The same calculation can be done in terms of Dyck paths. If we define hmin as the 
minimum height in the segment (*, j) 

hmin Min-f/ij, hij^\^ ■ ■ ■ 1 hj^i (4.14) 

we recover the generating function M(m; w): 

M{u]w) = ^;(^»-Vin) + (0-Vin)p(/^.^ fi.y (4.15) 

hi,hj 

There is an important conclusion coming from this exercise. The boundary two-point 
function of the bosonic field theory has a nonlocal lattice realization. Indeed, if we fix 
the segment, the end-points of the segment alone don’t fix the paths, one has to look at 
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Figure 7: The one-point function (14.Sh as a function of the lattice size L for the three values 
of w given in (14.Qh . From the hts on can read the values of A and A{w). 

the whole segment and hnd the value of hmin for each conhguration. In Sec. 5 we discuss 
the implications of the connection between the generic function M{u-,w) and the two-point 
function fl4.2p . 

The results described up to now can be generalized. Take again the segment {i,j) and 
denote u = j — i. Consider the generic function: 

M{u; wi, W 2 ) = Y. hj). (4.16) 

{hi} 

we conjecture that in the large L limit, we have 


M{u]Wi,W2) 


C{wi,W2) 

.yA(ji;i)+A(«)2) ’ 


(4.17) 


where C(l,l) = 1. Obviously we recover fl4.13p if wi = W 2 = w. Numerics suggest that 
C(wi,W 2 ) = A(wi)A(w 2 ), where A(w) is dehned by (14.81) . 

The conjecture fl4.17p was checked using Monte Carlo simulations for several pairs of 
values (wi, 102 ). In Figs. 8 and 9 we show the results for the pairs (l/\/3, y^ 2/3)) respectively 

(^2/3, y^ 3/4). We have taken several values of L: 300, 600, 1200, 2400, 4800 and 9600. 
One can see that the data are independent on L. Moreover, the u dependence is correctly 
given by (I4.17p when using the exponents (14.9p . 

We are not aware of any held theoretical derivation of the two-fugacity generating function 
so our results stay as a conjecture. 
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Figure 8: The two-point function M{u]Wi,W 2 ) defined in 04.171) as a function of u for the 
values wi = l/\/3 and W 2 = y^ 2/3. The predicted result follows from 04.9j) and is given in 
the hgure. The lattice sizes are L = 300, 600,1200, 2400,4800 and 9600. The fitting was 
done for the largest lattice and in the region of the dashed box shown in the hgure. 



Figure 9: Same as Fig. [8] for the two-point function one M{u; 101 , 102 ) dehned in 04.17P for 
the values wi = ^2/3 and 102 = \/3/2 
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Figure 10: Particle densitity correlator as a function of u Correlation of density of particles 
for the periodic RPM. The data are for the lattice sizes L = 1200, 2400, 9600,19200 and 
38400. The plot is in a log-log scale considering only the odd lattice distances. 

5 Correlation functions in the Raise and Peel model 
and conformal invariance 

In Sec. 2 we have given the predictions on conformal held theory for the correlators of 
stochastic models. These predictions were made for local boundary operators. We have 
described in Sec. 3 the Raise and Peel model which is conformal invariant but in the stochastic 
base, the Hamiltonian acts in a nonlocal way therefore the task is to identify local operators 
in the stochastic base. In Sec. 4 we have shown that the vertex operators of the bosonic 
boundary held theory are indeed local but that the lattice model from which they can be 
derived is nonlocal. In this section we are going to clarify the problem and make clearer the 
connection between the previous sections. We discuss separately the periodic system and 
the open one 

a) Periodic boundary conditions. 

We consider hrst the particles-vacancy stochastic basis. A natural local variable is the 
particle density. One expects (I2.5jl to be valid and this is indeed the case. In Fig. HU] we 
show the density-density correlator < ri{0)ri{u) > obtained from Monte Carlo simulations 
for diherent lattice sizes {ri{u) is the density of particles minus 1/2). There is no lattice 
dependence and one obtains an exponent A = 0.983 which within errors is equal to 1, one of 
the exponent given by fl2.2p . This exponent is the one one expects on dimensional grounds 
and will also be obtained in the open system. 

Another local observable is the current density in the stationary state. It was shown in 
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Figure 11: Current density correlator as a function of u for the periodic RPM. The data 
are for L = 9600 

[12] that its expression for large L is: 



(5.1) 


One would expect therefore on the two-current correlator < J(0) J(m) > an exponent 2A = 2. 
Monte Carlo simulations (see Fig. fTTl) give an exponent equal to 1.028 which is close to 1 
and not to 2. We have no explanation for this observation. 

There is no obvious connection between the bosonic boundary held theory and the density 
of particles as well the current density observables. 

We consider now the Dyck paths stochastic basis. A natural local observable is the density 
R{L) of contact points {hi = 0) (see Fig. 4). There is an exact result for this quantity. It 
is based on a conjecture by de Gier [T3| for the number of clusters for any L. A cluster is 
a segment between two consecutive contact points. This conjecture was verihed by Monte 
Carlo simulations. The expression of the density for large L is; 


R(L) = |V3r(5/6)/U/3 


(5.2) 


This gives an exponent A = 1/3, again given by (12.2p . The same result is obtained using 
fl4.3p and fl4.7p with tc = 0. 

The dehnition of the two-contact-point correlation has to be considered with care. If 
one uses the SOS prescription discussed in Sec. 4, one has to use Eq. i fl4.2p for tc = 0 and 
get from 04.31) and 04.7p A = 1/3. We have checked this prediction using Monte Carlo 
simulations. This implies that one has to take into account not only the conhgurations in 
which hi = hj = 0 (j — i = u) but also all the conhgurations in which hi = hj = h-min. Those 
are the NCA conhgurations. It turns out that if we take into account only the conhgurations 
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with contact points {hi = hj = 0) the two-point fnnction vanishes in the large L limit bnt it 
does so in a neat, nnexplained. way: 


The physical interpretation of fl5.3p is the following one. The conditional probability to have 
a contact point at the distance u from a given contact point (that happens with probability 
(decreases like 

Before closing this section we give some results which can be obtained using the generating 
functions M{w) (Eq. fl4.8l) ). M{u;w) (Eq. fl4.131) 1 and M{u]Wi,W 2 ) (Eq. fl4.17p l dehned in 
Sec. 4. One obtains in this way properties of the Dyck paths conhgurations in the large L 
limit. 

From the one-point vertex M{w) (4.8) one obtains 


< >= 



M{w) 


11)=! 


P{h = n) 


^ ^h,n 


1 

n\ dw^ 


M{z) 


w=0 


It is easy to derive the leading behavior: 


< h >= 


271 


In L, 



<e>-<h>‘^= 


2\/37r-9 


In L, 


and 


(5.4) 

(5.5) 


(5.6) 

(5.7) 


P{h = 0) = A{0)L 
P{h = 2) = 71(0) ( 


-1/3^ p{h=l) = A{0) 

27 (lnL)2 9 InL' 
327r2 L^A 


3\/3 In L 
4:71 L^A’’ 


with ^4(0) given by fl4.1ip . 

From M{u;wi,W 2 ) (Eq. (I4.16p ) we get the leading behavior: 


< hi 


hmin ^ 


d 

dwi 


M{u-,Wi,W2) 


Wi=W2 = l 


V3, 

= — mu. 
271 


(5.8) 


(5.9) 


Notice that in the large u limit < h — hmin > coincides with (5.6). This implies that the 
average number of arches crossing a link is equal to the average number of arches entering 
the segment of length u from one side or the other. The double of this quantity is called 
valence bond bipartite entanglement entropy (see e.g. in [T4l 13]). 

An interesting quantity is the average probability to have a height {hi — hmin) = m at 
the site i and a height {hj — hmin) = n at the site j. Using (I4.16p and (I4.17p one hnds the 
leading behavior: 

Fm,n(w) ~ I ln(n)”^+lM-2/3. (5.10) 
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b) Open boundary conditions. 


We consider first symmetric boundary conditions and a stochastic basis with particles 
and vacancies in an open segment of length L {L even). Let u be the distance from one of 
the boundaries and p the density of particles from which one has subtracted 1/2 (the average 
value). Using Monte Carlo simulations [15] we have found 


P 


C 


71 

L sm{7iu/L) ’ 


(5.11) 


where C* is a constant. Using fl2.12p we hnd A = 1 in agreement with the result obtained in 
the periodic boundary case. 

We next consider the density of contact points Rl{u). An exact result is known in this 
case urn- 

( TV 

RL{n)=cA , (5.12) 

\Lsm[7ru/L) J 

where Ci = — \/3r(— l/6)/(6(7r)®/®) = 0.7531... . Comparing with fl2.12p one hnds A = 1/3 
in agreement with what was obtained in the periodic case. 

The case of nonsymmetric boundaries was studied in Ref. [IS]. The one-boundary Tem- 
perley Lieb algebra was used adding a boundary generator to the Hamiltonian (13.ip and the 
function F{cos{7vu/L)) in Eq. fl2.14p was identihed for both the densities of particles and of 
contact points. 


6 Conclusions 

We think that in this paper we have clarihed the consequences of conformal invariance in 
stochastic processes. The main ingredient is the realization that if one looks at space depen¬ 
dent correlators one deals not with bulk operators but with Neumann boundary operators. 
Since there are no left and right movers in this case, the critical exponents are half of what 
one would expect for a periodic system. Space and time correlators imply that one has to 
consider also bulk operators in the vicinity of the boundary. One conclusion is that one 
shouldn’t expect the standard relativistic space-time dependence + R of the correlators. 

This picture was verihed in the Raise and Peel model which has a simple, albeit nonlocal 
dynamics which allows Monte Carlo simulations on large lattices. We have checked in this 
way the predictions of conformal invariance when analytical results were not available. 

Because of the nonlocal dynamics, the identihcation of local operators in the stochas¬ 
tic base for which conformal held theory makes predictions, is not an obvious task. This 
identihcation is dependent on which basis was chosen. 

In the particle-vacancy basis, we have checked that the particle density is a bona hde 
local operator with scaling dimension A = 1. The current density operator on the other hand 
has a correlator with an unexplained behavior. In the Dyck paths basis, one can use the 
mapping of the model onto the SOS model which in continuum space and time is described 
by a bosonic free held theory. Using boundary vertex operators with a known behavior, one 
can go backwards and identify the local operators in the lattice model. 


19 










Looking at Dyck paths, the density of contact (zero height) points look as a natural 
candidate for a local operator. Examining the behavior of the one-point function, it looks 
that this is the case. One hnds a scaling dimension A = 1/3. Looking however at the 
two-contact point correlator, we discover that it vanishes at large values of L. The L and 
u (the distance between the two contact points) dependence is simple, it involves again the 
exponent 1/3, but the expression of the correlator remains a puzzle. The way to solve the 
problem comes from the bosonic held theory. If one takes the segment between the two 
points where the correlator is computed, one sees that for each conhguration one has a 
minimum height hmin- If one subtracts from each height one gets new contact points 
for which h — hmin = 0. If one takes into account now all the conhgurations having contact 
points including also of those of the new kind at the end of the segment one hnds a non¬ 
vanishing correlator with the proper critical exponent. This is the NCA (non-crossing arches 
correlation). 

Let us observe that from all the possible scaling dimensions (2.3) available for a c = 0 
theory, we found only two of them: 1/3 and 1. 

As a byproduct of our research, we made the conjecture (4.12) on the average probability 
to have two diherent subtracted (hmin is taken away) heights at the end of a given segment. 
This conjecture was checked in Monte Carlo simulations. 

As the reader might have observed, we have not touched on the space and time correlators. 
To our knowledge, exact results are not known and computer simulations on very large 
lattices have large hnite-sizes corrections depending on the initial conditions. We plan to 
look in the future at this problem. 

Before closing the paper, we would like to observe that with the exception of several 
versions of the Raise and Peel model, no other models with conformal invariance and a 
simple dynamics are known. This is an invitation for more research in this held 
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